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What can be learned about dark energy evolution? 
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Abstract. We examine constraints obtained from SNIa surveys on a two parameter model of dark energy in which the equation 
of state w(z) = P{z)lp{z) undergoes a transition over a period significantly shorter than the Hubble time. We find that a transition 

between w 0.2 and w 1 (the first value being somewhat arbitrary) is allowed at redshifts as low as 0.1, despite the fact that 

data extend beyond z ~ 1. Surveys with the precision anticipated for space experiments should allow only slight improvement 
on this constraint, as a transition occurring at a redshift as low as ~ 0.17 could still remain undistinguishable from a standard 
cosmological constant. The addition of a prior on the matter density Q. m = 0.3 only modestly improves the constraints. Even 
deep space experiments would still fail to identify a rapid transition at a redshift above 0.5. These results illustrate that a Hubble 
diagram of distant SNIa alone will not reveal the actual nature of dark energy at a redshift above 0.2 and that only the local 
dynamics of the quintessence field can be infered from a SNIa Hubble diagram. Combinations, however, seem to be very 
efficient: we found that the combination of present day CMB data and SNIa already excludes a transition at redshifts below 0.8. 
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1. Introduction 

The nature of dark energy is one of the most puzzling 
mysteries of modern cosmology. It is now widely accepted 
that our u niverse is exper i encing a phase of accelerated 
expansion ([Peebles & Ratral 12003I) . The evidence was first 
found us ing a type la supernovae luminosity vs . redshift di- 
agram dRiess et al. 1998 : Perlmutter et al. 1999b . but a num- 
ber of other observations now support this conclusion. In 
particular, estimations of the matter density of the universe 
generally lead to a low value, while the CMB anisotropies 
point toward a spatially flat Universe dLineweaver et al. 1997 ; 
de B ernardis et al. 2000j). WMAP data also require the pres- 
ence of dark energy (ISpergel et al.1 120071) unless one con- 
siders an unexpectedly low value of the Hubble parameter 
H < 50 km s 1 Mpc" 1 dBlanchard et al.ll2003l:lHunt & Sarkat 
2007b . The possible correlation of the C MB fluctuation map 



with surveys of ext r agalactic objec ts dFosalba et al.l 12003 



Corasaniti et al .1 l2005UPogosianll2005l ) also provides direct ev- 
idence, although with a limited significance level (< 3<x), 
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for the existence of an unclustered dark energy component 
in the universe, the detected correlation being explainable 
through the integrated Sachs-Wolfe effect, i.e., a time varia- 
tion of the gravitational potential, which is achieved only if 
the (baryonic + cold) matter density parameter Q m signif- 
icantly differs from 1. Finally, the shape of the correlation 
function on scales up to lOO/i" 1 Mpc which has been recently 
measured accurately dEisenstein et al] 12005b in combination 
with the CM B data advoc ates for the presence of dark en- 
ergy (Blanchard et al . 2006) in the framework of general rela- 
tivity. However, the nature of this dark energy has been the sub- 
ject of numerous speculations. The simples t model, which was 
originally proposed (in another context) by lEinsteinl d 1 9 1 7 ). is 
a pure cosmological constant A, a term on the left hand side 
of Einstein's equations. However, a cosmological constant can 
also be regarded as the contribution of the vacuum to the right 
hand side of the equation with a specific equation of state, i.e., 
a component with negative pressure Pa related to the energy 
density pa by the relation Pa = -p\. Indeed, quantum field 
theory predicts that the lowest energy state of any mode con- 
tributes to a vacuum energy density that behaves exactly as a 
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cosmological constant (see, e.g. lBinetruv 2000). A number of 
problems arise with this possibility, in particular the so-called 
hierarchy problem: the expected contribution is usually enor- 
mous, naive calculation gives p vac ~ 10 76 GeV 4 , around 122 
orders of magnitude larger than the present critical density of 
the universe. However, there exist mechanisms, such as super- 
symmetry, which allow one to reduce considerably the vac- 
uum energy density, but since supersymmetry is broken at a 
scale larger than 100 GeV one is still plagued with an enor- 
mous vacuum energy density of the order of 10 s GeV 4 . The 
usual explanation is then to say that there exists a yet unknown 
mechanism which ensures that the contribution of the vacuum 
energy density is zero. One is therefore left to explain the na- 
ture of dark energy, which differs from a cosmological con- 
stant, avoiding the extreme fin e tuning require d to obtain the 
observed dark energy density dBinetruvl 120001) . Most models 
that have been proposed so far (quintesse nce models) therefore 



rely on the idea that some scalar field (ICaldwell et alJ 11998) 



behave today like a cosmological constant, exactly as an other 
scalar field did during inflation. The most remarkable feature 
of quintessence models is that both the scalar field pressure 
Pq and energy density pq evolve according to dynamical equa- 
tions. Consequently, the so-called equation of state parameter, 
wq = Pq/pq, varies with time between 1 and -1, as the field 
evolves along the potential. In some extreme (and possibly ill- 
defined) models, this parameter can even take any arbitrary 
value, for example if one allows the density pq to t ake nega- 
tive v alues or a change in the sign of the kinetic term (Caldwell 
20021) . Other models involvin g scalar tensor theo ries also allow 
for such transient behaviour (Eli zalde et alj|2004l) . The detec- 
tion of such a variation would therefore be of great importance 
for our understanding of dark energy. 

The aim of the present paper is to study models with a 
rapid transition of the equation of state and to illustrate that 
in this case, the Hubble diagram of SNIa provides surprisingly 
weak constraints compared to the case of a smooth transition. 
In Sec. |2j we recall a few basic aspects of simple quintessence 
models, and the motivation for a convenient parametrization of 
the equation of state parameter wq(z) allowing rapid transition. 
In Sec. [3] we describe the analysis we perform, and state our 
main results. In Sec. [4] we discuss the crucial issue of the im- 
pact of the epoch of observation on the parameter estimation. 
We draw the main conclusions of our work in Sec. [5] 

2. Dark energy parametrization 

Historica l quintessence models rely on the idea o f a tracking 
solution dRatra & Peebles! Il988t IWetterichl Il988l) . which in- 
volves a scalar field evolving in an inverse power law poten- 
tial, V(Q) oc Q a , the proportionality constant being tuned so 
as to obtain the desired value of the dark energy density pa- 
rameter Qq ~ 0.7 today. The main feature of these models is 
that the pressure to energy density ratio, wq, remains constant 
both in the radiation era and in the matter era (with different 
values during each epoch), and that it tends toward -1 once 
the quintessence energy density dominates. The value of the 
parameter wq depends on the power law index of its potential. 
From existing data, it seems that only values close to wq 1 



today, or even possibly lower, are acceptable dCaldwelll l2002; 
Melchio rri et al. 2003 ). Note however that in the latter case, 
values of wq below -1 cannot be obtained naturally through a 
standard scalar field. Single power law potentials suffer from 
the fact that, once quintessence dominates, the wq parameter 
approaches the asymptotic value -1 very slowly so that today, 
if the quintessence density parameter Qq is close to 0.7, then 
wq is still far from the value -1, contrary to what most analy- 
ses suggest. In order to avoid this problem, one has to add extra 
features in the potential, such as a rapid change in the slope of 
the potential o r a local minimum, su ch as in the SUGRA model 
proposed by iBrax & Ma rtin d!999h . Many other possibilities 
have been proposed since then (see for example references in 
Brax et al.ll2000l:IPeebles & Ratrall2003l) . 



On the other hand, without precise ideas about the correct 
quintessence model, it has become natural to adopt a more phe- 
nomenological approach in which one parametrizes the func- 
tional form of wq(z) which exhibits the main features described 
above. 

The simplest model of quintessence (in the sense that it 
introduces only one new parameter as compared to a ACDM 
model) is to assume a constant wq. However there is little mo- 
tivation for constant wq beyond the economical argument and 
it is increasingly recognized that evolving wq should be in- 
vestigated with a minimal number of priors. In the absence of 
well motivated theoretical considerations one is left with the 
empirical option to examine constraints on the analytical form 
for wq(z). Most investigations have been based on expressions 
with one or two parameters. However, such expressions often 
vary with time in a relatively slow way and that rapidly varying 
expressions have to be examined as well. In other words, if one 
considers the typical time scale: 



w 

W 



(1) 



constant w corresponds to tq » t H where t H = l/H is the 
Hubble time, a smoothly varying expression such as the inverse 
power law potential corresponds to tq ~ t H and a more rapidly 
varying w correspond to tq <k tn, such as in the SUGRA 
model. Our aim is primarily to investigate constraints on mod- 
els for which tq <sc f#. This lead us to use the following model 
which allows arbitrary rapid transitions and in which the dark 
energy wq parameter evolves as a function of the scale factor a 
according to 



1 1 I (a 
WqW = + Woo) - -(Wi -w 00 )tanh Tlog — 

2 2 \ \a t 



(2) 



The w parameter goes from w\ at early times to Woo at late times, 
the transition occurring at a t . The transition occurs at redshift 
Zt = l/a t - 1 (a negative value of which corresponds to a tran- 
sition in the future) and lasts of the order of F l Hubble times; 
T is therefore a parameter describing the speed of the transi- 
tion : high values (» 1) correspond to fast transitions, in the 
limit r = oo the t ransition is instantaneou s. This expression 
was proposed by iLinder & Hutererl ([2005). The quintessence 
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conservation equation D^T^ = can be integrated exactly to 3. Analysis 



give 



p Q (a) = p Q 



30+»v)/(fl/ flt ) r + ( fl /a t r 



qm) 1 +(i/ flt r 



(3) 



Aw = (wi - Woo), 



(4) 



where we have set 

Wav = ^(yVi + Woo), 

and where corresponds to the values of the quintessence en- 
ergy density at the present epoch. Therefore, the model can be 
implemented without muc h modification in existing cosmolog- 
ical codes such as CAMB jLewis et al.l2 000). Dark Energy per- 
turbati o ns are implemen ted accordingly to iRiazuelo & Uzan 
d2002h : iBrax et all d2000h . i.e. the initial spectrum of the dark 
energy perturbation follows an attractor mechanism, exactly as 
the unperturbed part of the DE fluid does, and the DE initial 
spectrum is lost. One modification concerns the sound speed 
defined as P/p which is equal to the w parameter when w is 
constant. When w is not constant, as is the case for our model, 
its exact value is known since the analytic forms of w and p 
are known. This model has four parameters. As we said ear- 
lier, simple quintessence model have fewer: constant w models 
correspond to Wi = w m , with undefined a t and F. With an in- 
verse power law potential, one has w m = -1, Wi = -2 /(a + 2), 
the epoch of transition is fixed approximately by the constraint 
^Qfe) ~ ^m(zt) ~ 0.5 and the duration of the transition is 
larger than the Hubble time (it depends on how steep the po- 
tential is, that is, on a). For the SUGRA potential the first 
two above constraints on w; and Woo remain, whereas the lat- 
ter are modified: the epoch of transition for wq no longer nec- 
essarily corresponds to the scalar field domination (but corre- 
sponds to the epoch where the field reaches a local minimum 
in its potential), and the transition duration is usually much 
shorter. Various other models have different predictions con- 
cerning these parameters, but the above parametrization is suf- 
ficiently general to encompass a large number of already pro- 
posed models. 

One of these parameters, W{, is not expected to be as rele- 
vant as the parameter a t : because at early times Oq(z) is neg- 
ligible compared to fi m , the quintessence field does not play a 
crucial role, at least with respect to supernovae and CMB data. 
One should impose a value of w\ slightly lower than 0, in or- 
der to ensure that at early times Qq <k Q m . The reason is that 
with W[ - and a low z transition Oq would be close to its 
present value at the recomb ination epoch (or ev en nucleosyn- 
thesis, see references in Peebles & Ratrall2003l) . At low red- 
shift, this would lead to a dramat ic suppression of th e cosmo- 
logical perturbation growth rate dPouspis et al. 2003 ). In addi- 
tion, we found that this introduces additional changes in the C\ 
curve at high / (i.e., other than changes due to the modification 
of the angular distance). For these reasons we fix w\ = -0.2. 
Putting a constraint on W&, is less desireable since it implicitly 
selects a limited class of models, which do not seem excluded 
by the data. We have chosen the value Woo = — 1, which seems 
in agreement with the present data, and we focus on the two 
remaining parameters, T and a t which describe the transition 
experienced by wq(z) between its early and late behaviour. 



We focus here on constraints that can be set in the transition 
parameters Zt and F, and we set w\ = -0.2 and Woo = — 1 as 
explained above. Note that a pure cosmological constant be- 
haviour is obtained by considering large Zt with a sufficiently 
small transition duration (so that it does not last long after Zt). 

3.1. Supernovae Hubble diagram 

The luminosity distance is one of the ma i n sou rces of con- 
straint on the nature of dark energy (lAstien (1200 lh ~). We there- 
fore first examine what kind of constraints the Supernovae 
Hubble diagram allows. The number of well observed SNIa 
has rapidly increased in recent years and a significant number 

of supernovae above redshift one have been detected. In the 

I 1 J — L 

following we us e the last compilati on from Davis et al. (2007) 
of recent SNIa dWood-Vasey et all d2007l) : iRiess et all d2007l) : 
AstieretalJ d2006lV ). 

We have examined constraints on our three parameters, a t 
I / ( I +z t ), r and Qq. We use the publicly available Mon te Carlo 
Markov Chains code, cosmomc (ILewis & Bridld d2002l) '). using 
the modified version of CAMB described above. 2-D contours 
shown in FiguresQ]and|4]encompass 68% and 95% confidence 
levels (CL). In each figure, the third parameter is marginalised 
over. 

The upper graph of Figure Q] presents the allowed regions 
in the plane: a t - 1/(1 + Zt) versus the inverse duration of 
the transition F. The lower graph gives the constraints in the 
a t -Q.Q plane. Constraints on possible transitions appear very 
weak: only sharp transitions at very low redshift (z t < 0.1 at 
the two sigma level) are firmly excluded. Surprisingly, the data 
suggest a transitio n at low redshift, a tendency that has bee n 



noticed elsewhere (Bassett et al. 2004; Corasaniti et al. 2004). 



However the significance level is low and a cosmological con- 
stant remains consistent with the data at the 2 sigma level. 

While rapid transitions (corresponding to large F) are very 
weakly constrained, better constraints are obtained when a 
strong prior is set on Q m : with Q m = 0.3 we found that tran- 
sitions are acceptable at redshifts greater than 0.25. This im- 
provement is due to the removing of degeneracy breaking (no- 
ticed in the parameter space of Figure [TJ but remains modest. 
This means that the Hubble diagram of distant SNIa alone is 
insufficient to determine the nature of the dark energy at high 
redshift. 

In Figures [2] we show the effect of a transition on the 
magnitude difference between a fiducial model and the empty 
universe, all other parameters being fixed with their fiducial 
values. The top figure clearly reveals that a transition from 
Wi = -0.2 to Woo - -1 occurring at even moderately low 
redshift makes very little differences to the observable quan- 
tity. It is therefore not surprising that the constraints that can 
be set from the present day SNIa Hubble diagram are not very 
tight. The bottom figure illustrates the effect of changing Wj: 
changing w\ from to -0.6 produces changes that are small 
and easy to understand as the model becomes degenerate with 
the ACDM model as Wi tends to -1. For this reason, in the 
following, we concentrate our analysis on a t . 
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Fig. 2. a) Residual Hubble diagram with respect to an empty 
universe for models with a transition at epoch a t = 1/(1 + Zt) = 
0.0, 0.3, 0.6, 0.7 & 0.8 compared with a binned version of the 
SNIa compilation of lDavis et al. (2007). Note that as explained 



Fig.l. Contour constraints on the transition at epoch a t = in the text low values of a t all give similar curves as the ACDM 

1 /( 1 + z t ), and the rate of the transition Y. The results are inde- model - b ) Same quantity for models with a transition at epoch 

pendent of F when it is sufficiently large. Bottom figure shows fl t = 0.565 and w; = 0, -0.2, -0.4, -0.6 
the results in the (£2q, a x ) plane. 



We have redone the above analysis on a t for a simulated 
survey with the precision and statistics expected from space 
experiments. We generated 2000 supernovae distributed in 16 
bins in redshift between 0.2 and 1.7 according to Table 1 of 
Kim et al. I d2004l) completed by 300 nearby supernovae. The 
number of supernovae per bin fluctuates according to a Poisson 
law. For a given bin, the magnitude of the supernovae is taken 
from a Gaussian distribution of the mean value given by the 
standard concordance ACDM model, and the sigma fixed to 0.2 



magnitude. The resulting magnitude of each bin is obtained by 
a fit of the distribution of magnitudes and the associated error is 
added in quadrature with a systematic error of 0.02 magnitude 
and an offset error of 0.01 magnitude for the intercalibration 
between the two sets of data. The constraints inferred from this 
simulated sample again reveal that the transition epoch a t is 
moderately constrained: transitions at redshift as low as 0.5 (2 
cr CL) are still acceptable when a rapid transition (F > 2) is 
assumed. 
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Fig. 3. Angular power spectrum of CMB fluctuations for mod- 
els presenting a transition at epoch at = 1/(1 + Zt) = 0, 0.3, 
0.6, 0.7 & 0.8 (all other parameters fixed) compared with the 
WMAP 3 TT spectrum. Apart from the large scale (low () 
power, the only difference between the spectra is a shift in the 
spectrum as a result of the modification of the angular distance 
to the last scattering surface. 



The situation is therefore paradoxical: although space sur- 
vey precision improves the constraints by pushing the accept- 
able redshift from 0.25 to 0.5 (for rapid transitions), this last 
result is modest as a significant fraction of the high precision 
data provided by the experiment extends up to redshift ~ 2. The 
reason for this apparent paradox is clarified in Sec. [4] 



3.2. Constraints from the CMB 

Given that SNIa Hubble diagram hardly suggests the presence 
of a transition in the dark energy content of the universe, it is 
interesting to examine whether such a possibility could be ac- 
ceptable using additional constraints. CMB is kno wn to provide 



such constraints on the quintessence scenario, eg dSpergel et al 



20031: iDouspis etalj 120031: lOdmanet all 12004 . We therefore 
examine CMB constraints on the type of models introduced 
above, although we leave to a future work a full investigation 
of the constraints that can be set on this type of model. We use 
the WMAP 3 data-set, as well as CBI, VSA and Boomerang 
data a t small scales, and a version of the CAMB cosmological 
code (iLewis et al.l l2000) that we have modified. Modifications 
of the code are straightforward since its public version includes 
models with constant wq in which we have implemented the 
energy density pq(z) and the pressure Pq(z) as a function of 
redshift. Our ansatz for the equation of state parameter Wq(z) 



0.8 - 




Fig. 4. Contour constraints on the transition epoch ajao = 
1/(1 + z t ) versus the value of the dark energy density param- 
eter at present time, Qq from the CMB alone (upper) using 
WMAP3 , CBI, VSA and Boomerang data. Constraints on the 
same quantities when combined with supernovae data (lower). 



allows us to integrate the conservation equation to obtain an 
analytical form forpQ(zfl 

The angular power spectrum of CMB fluctuations in the 
presence of dark energy is modi fied mainly throu gh the mod- 
ification of the angular distance dBlanchardll 1 9841) (see Fig. [3] 



1 One has also to make the distinction between the equation of state 
parameter wq = PqIpq and the "sound speed" squared c~ Q = Pq/Pq, 
which are identical when w is constant. 
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and Elgar0v & Multama ki (2007)). Although a strong depen- 
dence appears, this is partially lost through parameter degen- 
eracies which strongly weaken the final constraints. In addi- 
tion, ISW will contribute to lower levels as the transition is 
assumed at lower redshift, and this effect contributes to modify 
the angular power spectrum of the CMB fluctuations. We have 
investigated the CMB constraints on models with rapid transi- 
tions described by Eq. [2] F was set to 10. We have checked 
that varying F above 5 produces no appreciable differences. 
The other parameters that were left free are the baryon budget 
Da, the optical depth t, the Hubble constant Ho, the dark en- 
ergy density at the present day Qq, the index of the primordial 
spectrum n, the amplitude of fluctuations erg and the transition 
epoch a t = 1/(1 + Zt)' From the contours obtained in figure [4] 
one can see that the constraints that can be set on the transition 
redshift Zt from the CMB are rather stringent, Zt > 0.54 (lcr on 
one parameter) when Qq represents less than 50% of the total 
density. These constraints being slightly dependent on Hq we 
have also examined whether a combination of CMB and super- 
nova data allows to improve the transition epoch constraints, 
but although the SNIa data restricted the dark energy density 
much more around Qq ~ 0.7, the final constraints do not rep- 
resent a significant improvement: the final constraint shown in 
Figure H] is Zt > 0.66 (2<x on one parameter). Some different 
dark energy models could in principle lead to different conclu- 
sions in the case where the sound speed varies in a way that 
significantly affects the integrated Sachs-Wolfe effect on large 
angular scales, even though no such model was found in our 
analysis. Clearly, better constraints could be obtained from ad- 
ditional data of cosmological relevance, but this is beyond the 
scope of the present paper. 



between the two models, just as we are able to distinguish be- 
tween a ACDM with Qq = 0.5 and a flat Einstein-de Sitter 
model today. Now, are we able to distinguish today between 
these two models, which differ only in z > Zt? Surprisingly, 
the answer is no if one considers supernovae data only, as is 
convincingly illustrated by Fig|2] The explanation of this ap- 
parent paradox is as follows. Present data favour dark energy 
because high redshift supernovae are dimmer than expected in 
a flat Einstein-de Sitter universe. This is usually expressed as a 
difference of magnitude between the two models one considers 
for some standard candle at some redshift, the exact value of 
which depend on the quality of the data. The magnitude is es- 
sentially the logarithm of the luminosity distance as a function 
of the redshift. Let us define d£(z) and d^ ds (z) the luminosity 
distance as a function of the redshift in a ACDM model with 
Qa = Q m = 0.5 today, and in a flat Einstein-de Sitter model. 
Let us assume these two models can be distinguished. Let us 
now consider d^(z) and (ig(z) the luminosity distance vs. red- 
shift relation in a ACDM model with Q A = 0.7, Q m = 0.3 
today, and a dark energy model with Qg = 0.7, Q m = 0.3 to- 
day, with wq experiencing a sudden transition from to - 1 at 
z - Zt- An observer at z = Zt would therefore measure either 
d^(z') or df ds (z'). The epoch corresponding to a redshift of z' 
measured by an observer at z t corresponds to a redshift z given 

by 



z = (1 + z t )(l +z') - 1 = Zt + z" + z t z', 



(5) 



measured by an observer today. Let us define d t as the luminos- 
ity distance of the observer at z — Zt as seen from today. The 
exact value of d t does not matter here, but it can be computed 



4. Observing the Universe at z = and at z = 0.3 

In order to distinguish quintessence models from a pure cos- 
mological constant, it is crucial to be able to track the dark 
energy evolution as early as possible. The main impact of dark 
energy comes from its influence on the expansion rate of the 
universe. An important question is therefore until what epoch 
the dark energy density plays a role in observable quantities, 
and as a corollary, until what epoch one can hope to recon- 
struct either its energy density or its equation of state param- 
eter. As we have seen in section [3] the SNIa Hubble diagram 
poorly constrains a possible transition epoch in the equation of 
state of the dark energy component. As we have stated, this ap- 
pears somewhat paradoxical as data extending up to redshift 
2 fail to reveal a transition occurring at redshifts as low as 
0.25, at which the dark energy component is still dominant. 
Indeed, in a model with Qq = 0.7, Q m = 0.3 today, with 
wq constant and equal to -1, the matter to dark energy tran- 
sition, defined when Qq(z) = Q m (z) = 0.5, occurs at redshift 
z e = (QQ/Q m )5 - 1 ~ 0.33. Let us now consider two alter- 
natives. First, we can consider a pure cosmological constant 
model, with wq = -1 also at early times. Second, we can con- 
sider a model where wq ~ for z > z t = z e - In the first case, 
one has a usual ACDM model, whereas in the second case, one 
has a model close to a flat Einstein-de Sitter model at epoch 
z > z e - An observer at z = Zt should easily be able to distinguish 



1 



1+Zt 



- 1 



(6) 



Luminosity distances do not add but are proportional to comov- 
ing distances. The luminosity distance at redshifts above z t is 
therefore given by 



diiz) = l^d x + - 1+ , Z N ^(z'(z)). 



1 + z t 1 + z'{z) 
1 +z , 1 +z 



d L& = + : : ',;x d V(z)) 

L 1 + Zt 1 + z'iz) 



z'iz) 



(7) 
(8) 



in the two models. Now, in term of difference of magnitude 
between the two models, this is translated into 



Am(z) oc log 



d > + T&k d L&W 



whereas for the observer at z,, this simply gives 



Am(z') oc log 



dtiz') 
df dS iz')' 



(9) 



(10) 



Two differences arise here. First, the redshift range which to- 
day's observer must use in order to distinguish between the 
two models is larger than that of the observer at z t - If the lat- 
ter must collect data between z' = and z' = z», say, then 
the former must observe supernovae between z(0) = z t and 
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z(z«) = z t + Z* + ZtZ* > Zt + Z*. Second, the magnitude dif- 
ference between the two models is strongly attenuated for the 
observer today because of the presence of the extra term d x both 
in the numerator and the denominator of Eq. ©. This behavior 
is somewhat surprising at first but can be understood from the 
expression of the comoving coordinate r : 



r(z) 



dz 



)o H(z) 

where H(z) follow the Fridmann-Lemaitre equation: 
H\z) = %7iG Pm +p Q 



(ID 



(12) 



(for k = 0). Let us assume that a rapid transition occurs at z = z t 
between w = -1 and w ~ 0; it is straightforward to compute 
the relation between the Hubble constant H(z) in such a model 
to the Hubble constant H(z) in the standard model: 



H\z) 
H 2 (z) 



1 + 



l.-flo 

no(i+z0 3 



1 + 



l.-Qo 

n„(i+ Z ) 3 



(13) 



It is clear from this expression that the values of the Hubble 
constant in the two models do not differ by much when the 
transition redshift is not close to zero. As the integrand in the 
calculation of r is higher at low redshift, the limited difference 
in H translates to a small difference in r. For instance a tran- 
sition at a t ~ 0.7 results in a 5% decrease in r at redshift 1 
corresponding to Am ~ 0.1 as can be seen in figure [2] 

The net result is that while both models are easy to distin- 
guish at z = 0.3, this is no longer the case at z = as seen in 
Fig.H 

5. Conclusion 

We have investigated a class of models that undergo a rapid 
transition in the equation of state of their dark energy com- 
ponent. In order to establish the constraints that can be ob- 
tained on the characteristics of the transition we have focused 
our study on a class of models in which dark energy transits 

rapidly between w ~ and w 1 . We found that the duration 

of the transition cannot be constrained when it is shorter than 
the Hubble time. More surprisingly we found that SNIa Hubble 
diagram does not constrain this type of scenario very much, as 
even with the data expected from space experiments a transition 
can still be allowed at epochs when the dark energy density rep- 
resents 40% of the density of the Universe. This suggests that 
the SNIa diagram is poorly sensitive to dynamics of dark en- 
ergy at redshifts above 0.5. On the contrary, we found that ex- 
isting CMB data, in combination with SNIa data, provide tight 
constraints on this type of scenario, allowing rapid transitions 
to happen only at redshifts beyond 0.66, when dark energy rep- 
resents less than 10% of the total density of the Universe. This 
illustrates the importance of combinations of various data in 
order to accurately constrain the evolution of dark energy. 

Acknowledgements. The authors acknowledge Programme National 
de Cosmologie for financial support during the preparation on this 
work. The authors also acknowledge N.Aghanim for her comments. 
M.D. would like to acknowledge the French Space agency (CNES) 




0.2 



0.4 



0.6 



^A /0 mat Z obs 



0.8 1 1.2 

= 0.0, W; = 0. 



1.4 

£1 




£l./£l ,z. =0.3, w. = 0. 

A mat obs ' 1 



Fig. 5. Comparison of Q m -Qa fit for the two models discussed 
in Sec. [4] for an observer at z = (after the transition in the 
quintessence model) and z = 0.3 (during the transition). At 
Z = 0.3 the two models are very easily distinguished as one of 
them corresponds to a pure matter model, whereas one corre- 
sponds to a cosmological constant model with £2a = 0.5. On 
the contrary, at z = both models look like a cosmological 
constant model, and cannot be distinguished at the 2<x level. 
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